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ABSTRACT 

Observations of the intra-cluster medium (ICM) in galactic globular clusters (GCs) show a 
systematic deficiency in ICM mass as compared to that expected from accumulation of stellar 
winds in the time available between galactic plane crossings. In this paper, we reexamine the 
original hypothesis of Scott and Durisen that hydrogen-rich explosions on accreting white dwarfs, 
classical novae (CNe), will sweep out the ICM from the cluster more frequently than galactic 
plane crossings. From the CNe rate and stellar mass-loss rate, this clearing mechanism predicts 
that « 0.03 Mq should be present in < 10 5 Af Q GCs. We model the expanding remnant made 
from the 1O~ 4 M nova ejecta and show that it escapes long before it has cooled. We discuss 
the few positive ICM measurements and use a Monte-Carlo simulation of the accumulation and 
CNe recurrence times to reveal the possible variance in the ICM masses for the higher mass 
(> 5 x 10 5 Mq) GCs. We find that nova shells are effective at clearing the ICM in low-mass 
GCs (< 10 5 Mq), whereas higher-mass clusters may experience a quiescent time between novae 
long enough to prevent the next nova shell from escaping. The nova clearing mechanism will also 
operate in ultra-faint Milky Way satellites, where many upper limits on gas masses are available. 



Subject headings: galaxies: dwarf 
waves — stars: mass-loss 



globular clusters: general — novae, cataclysmic variables — shock 



1. Introduction 

The stars in a galactic globular cluster (GC) 
lose mass during their evolution, injecting ma- 
terial into the GC that accumulates over time 
to make an intra-cluster medium (ICM) of mass 
Micm- The most robust mechanism for clear- 
ing this matter is ram pressure stripping during 
the GCs passage through the galactic disk every 
~ 10 8 - 10 9 yr, implyin g M ICM ~ 100 - 1OOOM 
(Tayler & Wood 1975). However, searches for 



the ICM have yielded upper limits (Birkinshaw 
et al. 1983 Roberts 19881) or detections at levels 



more than 10 times smaller than implied by this 



mechanism ( 


Smith et al. 


1976 


Freire et al. 


2001 


van Loon et al. 


2006). Recent searches for dust 



in GCs confirm this paucity (Lynch & Rossano 
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1990l|Knapp et al.|1995[|Origlia et al.|1996[|Penny 



et al |1997[|Hopwood et al.|1999[|Matsunaga et al 
2008). Spitzer observations place upper bounds 
on dust masses 10 — 100 x below that expected 
), and, where detections are 
20061 120081), the values are 



(Barmby et al. 2009 



made (Boyer et al. 



much lower than expected. 

It's unclear if the diversity in GC escape ve- 
locities and the possibility for fast stellar winds 



(Dupree et al. 20091 can explain all of these ob- 
servations by simply having the material directly 
leave the cluster. A more robust solution for all 
GCs would be a more frequent ICM clearing mech- 
anism than disk passages. Though many mech- 
anisms (e.g. pulsar winds (Spergel 1991), stel- 



lar collisions (Umbreit et al. 



20081, ram pressure 



stripping by galactic halo medium ( Priestley et al 



2010)) are discussed, one appears to be robust 
(but often neglected) : ICM clearing from classical 
novae (CNe) explosions (Scott & Durisen 19781. 



Scott & Durisen (1978) showed that CNe had ex- 



1 



plosion energies, E nova , and recurrence rates ade- 
quate to unbind the accumulated ICM, but the 
astrophysical uncertainties (at that time) made 
definitive statements about this mechanism diffi- 
cult. We show here that the improved knowledge 
of mass loss from the stellar population; coupled 
with better constraints on the CNe rates and en- 
ergetics reveals that CNe will clear the ICM and 
do so with a frequency that explains the ICM ob- 
servations. 

The mass loss rate from the old stellar popula- 
tion in a GC, Micm, is now observationally well 
constrained. The main sequence turnoff mass is 
w 0.78 - O.88M potter et al.||2009| iThompson 



et al. 2010) in an old (> 10 Gyr) GC, and white 



dwarfs (WDs) are born with an average mass of 
w 0.53 ± 0.01M© ( |Kalirai et al||2009| ) implying 
« O.3M of material placed into the ICM^We es- 
timate the WD birthrate from the calculated spe- 
cific evolutionary flux for a Salpeter initial mass 
function (IMF) of 1.7x 10" 11 L^ 1 yr" 1 (from |Buz- 
zoni et al. (2006), where a steeper of shallower 
IMF only changes this number by ±20%) giving 
Micm ~ 5 x 10~ 12 Mq yr _1 Lq 1 . Assuming that 
each CNe clears the accumulated ICM, the typi- 
cal ICM mass is Micm ~ 0.03 M Q for our fiducial 
(see §2) CNe rate of 20 yr" 1 per 10 11 Lq. The 
ICM binding energy in a GC of mass Mqc and 
half-mass radius Rh is 



10 45 



ergs 



M, 



GC 



5 x 10 5 M, 



© 



M 



ICM 



1 pc 

O.O3M y \~R^ 
(1) 

similar to the CNe energy release, E nova , as first 



highlighted by Scott fc Durisen (1978) 



These arguments motivate our re-consideration 
of the dynamics of CNe shell evolution in §3, where 
we show that the nova shells escape the cluster be- 
fore suffering significant cooling. In §4, we apply 
this model to all galactic GCs, and discuss the 
possibility that the current ICM knowledge can 
constrain the underlying population of binaries 
that cause CNe, the cataclysmic variables (CVs). 
We conclude in §5 and mention the impact that 
CNe will have on the ICM in the dwarf spheroidal 
galaxies of the Milky Way's halo. 

1 We do not distinguish between the matter lost on the way 
to the Horizontal Branch (see Dotter 2009 for a recent dis- 
cussion) versus that lost on the asymptotic giant branch. 



2. Classical Novae properties and rates 

CNe result from the unstable thermonuclear 
ignition of accumulated hydrogen on a WD in 



a mass-transferring binary ( |Gehrz et al. 1998) 



The rapid nuclear burning drives the ejection of 
most of the accumulated material, leading to out- 



flows of M„ 



10" 



io- 



Mq at velocities of 



Kj ~ 1000 km s" 1 (see |Starrfield et al.| ([1972]; 
Prialnik ( |1986 )). The typical kinetic energy is best 
inferred from observations late in the CNe evolu- 



tion (typically > 10 years later (Cohen & Rosen 



thal|1983||Cohen|1985||DownesfeDuerbeck|2000| )j 

where V e j ~ 700 — 1200 km s _i is measured for the 
shell, motivating a fiducial CNe explosion energy 
of 



-S'nnva. — 



= 10 45 erg 



10- 4 M, 



o 



1000 km s" 1 
(2) 

which we use throughout this paper. The calcu- 
lated accumulated masses depend on the mass ac- 
cretion rate, WD core temperature and accreted 



metallicity (Prialnik & Kovetz (1995); Townsley 



fe Bildsten| ( [2005| ); |Yaron et al.| ([2005]) and, for 
a 'typical' case yield M ej « 3(10) x 10" 5 M for 
Mwd = 1-0(0.6) Mq, consistent with the observed 
range (Gehrz et al. 19981. This variation in CNe 



ejecta masses motivates our later consideration of 
the impact of a range of E nova . 

The CNe rate in GCs is not known, as only two 
CNe have been observed in galactic GCs (see his- 



torical discussions in Shara et al. (2004); Shafter 



& Quimby (2007)), and three in extragalactic 
GCs (one of M87 (iShara et al.||2004|), and two 



of M31 (Shafter fc Quimby 



2009 Pietsch et al. 20101) 



2007 Hcnzc et al. 



Shafter & Quimby 



(2007) discussed this challenge, suggesting that 
the CNe rate per unit stellar mass in GCs may 
only be slightly higher than in the field, but is 
likely not much lower. Another way to make an 
estimate is to take advantage of the explicit con- 
nection between CNe and the observed CV pop- 



ulation. Townsley & Bildsten ( 2005 ) showed that 



the majority of galactic CNe arose from CVs with 
-Porb < 8 hours undergoing the conventional mass 
transfer scenario of disrupted magnetic braking 
pOOll ) 



see 



Howell et al. 



distant galaxies, 20 yr 



The meas ured CNe in 
per 10 11 Lq (jGiith et al 



2010 ) then implied a CV population per unit mass 



comparable to that seen in our local galactic disk. 
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The combination of HST and Chandra obser- 
vations of GCs are revealing their CV populations 



(Edmonds et al. 


2003; Dieball et al.||2005 2007 


Bassa et al.|2008 


Dieball et al.|2010[ |Huang et al. 


2010 
their 
appe 
diall 


), allowing for much better raw sta 
incidence. However, these CV po 
ar to be a combination of those fror 
binaries and dynamical interactions ( 


tistics on 
pulations 
q primor- 
Pooley & 



Hut 2006 Stefano & Rappaport 1994 ), and may, in 



the end, depend on the GC central stellar density, 
po ■ Haggard et al. ( 2009 ) discuss the situation in 



u> Cen, where they conclude a lower CV incidence 
per unit mass than in the field. Given this current 
diversity of answers at present, we simply adopt 
the field C Ne rate of 20 yr" 1 per 10 11 L© ( [Girth 
et al.||2010[ ) for our workQ 



3. Remnant Evolution and Escape 

We now show that the kinetic energy injected 
by a CNe, E nova , is adequate to unbind the ICM 
accumulated since the prior explosion. Although 
local ICM structure will influence the shell evolu- 
tion, our aim here is to show the ease of unbinding 
the ICM. We define a characteristic (and uniform) 
ICM number density n = iMiQyi/m p AiTR\ ~ 
0.1 cm~ 3 , an overestimate since some ICM mass 
will be beyond Rh- We treat the evolution using 



a three-stage model (Spitzer 19781, starting with 



the brief free-expansion phase that lasts the time 

5/6 , p N -1/2 



tj = 120 yr 



10- 4 Mt 



V 10 45 ergs 



vCm 
(3) 

needed to sweep up the M e j of ICM, creating a 
reverse shock. At this point the evolution transi- 
tions into the Sedov-Taylor phase and ends when 
the radiative cooling of the post shock material 
becomes significant. 

The simple case of a strong explosion in a ho- 
mogeneous medium is given by the Sedov-Taylor 



solution ( Zel'dovich & Raizer 2002 1 . The radius of 



the shock front, i? s , velocity of the shock front, V s , 
and temperature immediately behind the shock, 



2 We can be confident that a GC will contain at least one CV 
since mass-transferring CVs have a birthrate of one every 
400 years per 10 11 M Q {Townsley fc Bildsten|2005[ |. If we 
assume an active lifetime of 3 Gyrs, then one active CV 
requires M GC > 10 4 Mq . 



T s are respectively, 

( 



R s (t) « 2.3 pc 



V s (t) « 92 km/s 



V 10 45 ergs 



1/5 



n 



-1/5 



t 



T s (t) = 1.1x10 s K 



V 10 45 ergs 



1/5 



Vcm 6 J 



10 4 yr 
(4) 

-V5 ( t 



\ 10 45 ergs 



2/5 



10 4 yr 

(5) 



Vcnr 



-2/5 



10 4 



(6) 

where t is the time since explosion. These numbers 
show that the shell will be moving at a speed in ex- 
cess of the GC escape velocity, V csc ~ 30 km s _1 , 
even 10 4 years after the explosion, consistent with 
our energetic arguments. The time it takes the 
shell to escape the cluster, t osc , is simply the time 
for the shell to reach Rh with V s > V csc . These are 
shown by the solid lines in Fi gure [Tj where we have 



used the Plummer model (Binney & Tremaine 



2008) to estimate V csc . 

The only remaining question is whether the 
remnant radiates and cools prior to escape. To 
answer that, we calculate the cooling of the rem- 
nant using collisional ionization equilibrium cool- 
ing rates for a range of metallicities from |Gnat fc| 



Sternberg ( 2006 ) . The mass of the shell is heavily 
concentrated near the shock front, therefore we ap- 
proximate the remnant as being spatially isother- 
mal at the temperature just behind the shock, 



P4)' 

Vcm J / 



V3T s (i). From Cox (1972), the energy loss rate is 



dE 
~dt 



(T,n) 



A(T)n 2 (R)A7rR 2 dR 1 (7) 



where A is the cooling function (in erg cm s ), 
so An 2 is the energy loss rate per unit volume. We 



adopt Cox ( 1972 ) 's compaction parameter 

'n(R) 



A = 



-ttR: 



no 



4irR 2 dR, (8) 



which has the value A ~ 2.3 for a strong shock. 
This yields 



dE 

~dt 



= -Ann 



; wRi(t) A 



(9) 



an equation that gives the energy as a function 
of time, defining i coo i, the time it takes for half 
the energy to be lost. We also compute the longer 



non-equilibrium cooling times for the case where 
the gas cools on a timescale shorter than the re- 
combination times of the various ions, causing a 
relative over-ionization as compared to the equi- 
librium case. We therefore distinguish cooling 
times that are calculated from the equilibrium 
cooling efficiencies, i C ooi,cq, and two types of non- 
equilibrium cooling efficiencies - isobaric, i coo i.ib 5 
and isochoric, t coo i ic . The upper curves in Figure 
[T] show the different cooling times for a range of 
R h for M G c = 5 x 1 O 5 M (and ICM particle den- 
sities ranging from no = 0.06 — 4 cm~ 3 implied 
by a fixed Micm = 0.03M Q ). The isobaric and 
isochoric cooling times are similar so we only plot 
^cooi.ic- All of these times are much longer than 
the escape time, making it clear that the shell will 
escape in the Sedov phase and radiative cooling 
will be an important energy loss mechanism only 
after the shell leaves the GC. 

4. Application to Galactic Clusters 

We now apply our model of ICM accumulation 
to the Milky Way GCs dHarris||1996] ) . We assume 
each GC has accumulated Micm ~ 0.03 M Q from 
stellar mass loss since the last clearing event, and 
that the characteristic particle density, no, is uni- 
form. For each GC, we compute the equilibrium 
cooling time, the shell radius at the cooling time, 
-Rcooi, and the shell's velocity, Vh, when it was at 
the half-mass radius (i.e. Vh = V s when R s = Rh)- 
We plot these values in Figure [2j showing that 
Vh > V csc in all cases, and that cooling occurs 
when R s 3> Rh- Escape is the outcome in this 
simple formulation. 

Since we don't expect the CNe to be strictly pe- 
riodic, we performed a Monte Carlo simulation for 
10 s years allowing for variation in CNe recurrence 
times and drawing the nova energies from a log- 
uniform distribution ranging from 2.5 — 10 x 10 44 
ergs. For each GC, we calculate the expected 
number of CVs (assuming 1 CV per 10 4 M Q , see 
footnote [2]), and generate a nova explosion history 
assuming an average CNe recurrence time in the 
range (3.3 — 6.7) x 10 5 yrs, centered on the recur- 
rence time calculated from the CNe rate, 5 x 10 5 
yrs. 

We calculate how large an effect these variances 
will have by calculating how often a GC is in a 
state where Micm > 0.03 M . During the 10 s 



FT 



I I I I I I I I I I I I I I I I I I I I H 



^cool,ic 



s * ^cool,eq 




1.0 1.5 

Rh (pc) 



2.0 



Fig. 1. — Timescales as a function of half-mass ra- 
dius for a GC of mass 5 x 10 5 M©. The top three 
lines are various cooling times. The topmost dot- 
ted line is the non-equilibrium cooling time from 
isochoric cooling rates. The next two dashed lines 
are the equilibrium cooling times. The heavy line 
is corresponds to our fiducial Micm = 0.03 M Q 
and the lighter line to a lOx higher Micm for com- 
parison. The lowest two lines are the escape times 
for these same ICM masses, again with the fiducial 
in bold. 



years, there is a 12% chance that the GC will have 
twice the fiducial Micm, a 4% chance for 3x, and 
a 0.4% chance to have 5x the fiducial ICM built 
up - a prediction nearly independent of GC param- 
eters. Shell escape can become difficult for some 
galactic GCs if Micm is increased by a factor of 
10. Figure [2] looks qualitatively the same in that 
case, but the line representing Vh — V csc would 
now cross between the triangles on the y-axis and 
the x-axis range would shrink to 6 — 14. 

The maximum quiescent time was ~ 10 x the 
average time between novae, allowing for a very 
large buildup of Micm, potentially to values where 
Eh > E nova , forcing us to consider the possibil- 
ity that some GCs may, due to the randomness 
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Fig. 2. — Outcomes for CNe driven remnants in 
Milky Way GCs. The x-axis is the shell's radius at 
the cooling time, i? C ooi, divided by the GCs half- 
mass radius, R^. The y-axis is the ratio of velocity 
of the shell at Rh to the GCs escape velocity, V csc - 
The shapes indicate GCs with ICM observations: 
47 Tuc is the up triangle; M15 the square; and M2 
the down triangle. The starred cluster is M80, a 
Milky Way GC with an observed CNe (T Sco in 
1860). 



of the process, experience a "runaway" accumula- 
tion. We determine the probability for this out- 
come by computing the average time it takes for 
this condition to be reached and assuming that the 
GC has uniform probability of being at any point 
in its 10 s yr orbital period (longer orbital periods 
will increase the chance of seeing a runaway) . For 
Pa = 10 3 L Q /pc 3 and Rh = 2.0pc, we find that 
GCs with Mqc = (1.5-5) x 10 5 M Q have an even 
chance of experiencing a runaway, where the ex- 
act value depends on the nova energy distribution. 
This probability increases very quickly with mass, 
so more massive clusters are much more likely to 
undergo a runaway while less massive GCs will 
almost never have one. This probability is also 
sensitive to po, with more collapsed clusters much 



more likely to undergo runaways. Although we 
refer to it as a "runaway" , it may not be unmit- 
igated mass accumulation, as the asymmetries of 
the ICM will allow the CN to unbind some mate- 
rial. However, the stochastic nature of this process 
inhibits our ability to make a secure prediction for 
a specific GC. 

There are also some instances (e.g. M15) where 
the shell escape time, i osc , is so long that it does 
not escape before another CN occurs. The out- 
comes from these colliding remnants are more 
complex, cither leading to a net heating of the 



ICM that creates a wind like solution (Scott & 
Durisen||1978|) or more intense cooling that keeps 



the material bound. In either case, novae clearing 
is most effective in the smaller, less dense clusters, 
where the shell safely escapes before another nova 
goes off and the GC is in no danger of a runaway 
scenario. 

There have been few positive ICM measure- 



ments in GCs. Barmby et al. (2009) derive Micm 



upper limits of 0.04 — 1.4 M Q from dust mass up- 
per limits. These upper limits are (with the excep- 
tion of M2) consistent with observed Micm values 
fro m|Freire et al.| pOOlj), |van Loon et al.| ( |2006[ ) , 
and iGrindlay & Lillerl (I1977I). We collect the four 



reported Micm measurements in Table [T] along 
with other useful GC parameters and timescales 
we derive. The ranges in t osc and t coo \ come from 
taking our fiducial Micm = 0.03 M in the clus- 
ters for one bound (the upper bound on i coo i and 
the lower bound on i esc ) and the measured Micm 
for the opposite bounds. Selection effects will bias 
likely detections towards GCs with higher Micm 
values. According to our calculations, the high 
central densities and masses of these clusters place 
them in the runaway regime. 

We have assumed that the mass lost from stel- 
lar evolution is a smooth process at the rate of 



Mi 



ICM 



5 x 10" 



M 



,•-1 



l:. 



However, 



stochastic effects in the mass loss should be ex- 
amined as well. We start by considering the most 
extreme case, where the mass is lost in a sudden 
event (e.g. He core flash, thermal pulse on the 
AGB). Since the WD birth rate is one-tenth the 
CNe occurrence rate, such an event would occur on 
an average timescale of 10i rcc . Though certainly 
the case that the CNe immediately following this 
event would not be adequate to unbind it (both 
because of energetic arguments and relative filling 
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fractions), such a stellar ejection "chunk" would 
be shocked, on average, by w 10 CNe. Whether 
this proves sufficient to unbind it depends on the 
location and time to cool between shocks; which 
is beyond the scope of our work here. 

An additional possibility is wind-driven mass 
loss that is so strongly dependent on the stellar 
luminosity (e.g. M oc L a ), that it makes Mtcm 
time dependent on the i rec timescale. As an ex- 
ample, consider mass-loss near the tip (Xtrgb) 
of the RGB. It takes 3 x 10 6 yrs to evolve from 
L = Ltrgb/2 to L = Ltrgb (Salaris & Cassisi 



1997), implying ss 30 stars in these final stages for 
a 5 x 10 5 Mq cluster. Let's consider the amount of 
mass one of these stars could lose during the final 
t rcc of its evolution to the expected total mass loss 
along the RGB, ~ O.2M ( jOrigliaet al.|2007b. For 
typical values of a ~ 1 — 3 ( McDona ld et al.|2 009 



Meszaros et al. 20091, this ratio is less than 10 



To get a substantial mass loss in the last t rcc re- 
quires a > 100, basically becoming equivalent to 
the sudden event model already discussed. So, we 
do not expect any of the current wind models or 
observations to create a time dependent Micm on 
the CNe recurrence timescale. 

5. Conclusions and other applications 

We have reconsidered the IScott fc Durisenl 



(1978) model of using CNe explosions to clear 



the accumulated ICM in galactic GCs, allowing 
for a more secure analysis of its efficacy. Using 
estimated CNe rates and the known stellar mass 
loss rates, we predict that novae clearing is highly 
effective in the less massive and less centrally con- 
densed galactic GCs (those unlikely to undergo 
a runaway), leading to a typical ICM buildup of 
« O.O3M that is consistent with observational 
upper limits. For those few GCs with positive 
ICM measurements, the Micm values are higher 
than our predictions, though still 1-2 orders of 
magnitude less than the ~ 300 M predicted from 
galactic plane crossing. We discuss in §4 whether 
this could simply be the result of the stochastic 
nature of CNe explosions and the large number of 
GCs observed, or from a systematic dependence 
on Mqc or p . In either case, it is clear that CNe 
have a large deleterious impact on the ability for 
the GC to build up it's ICM. 

A very robust mechanism to clear the ICM are 



Type la supernovae. The simplest assumption 
( Pfahl et al.|2009 ) is to use the rate expected from 
the mass, (5.3 ± 1.1) x 10~ 14 yr" 1 M" 1 ( |Sulli 
van et al.|2006 ) , yielding an average time between 
events of 2 x 10 7 yrs in a 10 6 M GC, shorter 
than a galactic plane crossing. If CNe are ineffec- 
tive in the intervening time, then the ICM mass 
would reach Micm ~ 100 M , but with a binding- 
energy orders of magnitude less than the 10 51 ergs 
of the SNe. Such outcomes may be more prevalent 
in massive GCs around elliptical galaxies, since 
there is no plane crossing to clear out accumulated 
ICM. When Micm ~ 100 M the explosion would 
occur in a density of n 100 cm~ 3 , much larger 
than expected in an elliptical galaxy, s» 10~ 2 cm~ 3 



( Stewart et al.|[l984 ). Application of our work in 
§2 finds that the velocity of such a remnant would 
be V s ~ 130 km s _1 at a time of 2700 years, when 
it reaches Rh = 1 pc. The cooling time is orders 
of magnitude longer, and so this shell will easily 
escape, and be quite bright (due to the high den- 
sity). Less than a few % of extragalactic GCs have 
shown emission lines consistent with these veloci- 
ties ( Chomiuk et al.|2008 |, raising the question as 
to whether these could be the first indicators of 
the expected la explosions in GCs. If such a rem- 
nant is detectable for 10 4 years, then we would 
expect to see one such case amongst all the GCs 
in a large elliptical galaxy like M87. 

Novae clearing can also be applied to the newly 
discovered ultra-faint (only rj 1O 4 M in stars, 
potentially so low that only a few CVs may be 
present, see footno te [2| dwarf sph eroidals (Be 



lokurov et al. 2007 Simon & Geha 2007), where 



ICM deficits have been reported ( Bailin & Ford 
2007[ )p] Just as in GCs, these dwarf spheroidals 
(dSph) have a mass-losing stellar population. 
However, the lower escape velocities (2 — 5 km/s 
from the s» 1O 6 M of dark matter in a half-mass 
radius of 100 pc) may allow the stellar winds to di- 
rectly leave the dSph. In cases where the mass ac- 
cumulates, the ICM binding energy is < 10 45 ergs 
for Mjcm = 0.03 Mq, so that a single CNe can be 



3 The large HI mass (co mparable to the stellar m ass of 
2 X 1O 5 M ) measure d by |Ryan- Weber et al.| | |2008| in the 
distant dSpH Leo T ( |Irwin et al.||2007[ l exceeds that pos- 
sible from stellar wind accumulation over 10 Gyrs, and is 
likely from continued infall of material. Such an inflow may 
also be the source of the 10% o f stellar mass in young (rj 
Gyr) stars ( de Jong et al.|2008 l. 
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Table 1 



GCS WITH ICM DETECTIONS 



Cluster 


M GC 


Rh 


[Fe/H] 


MiCM 


-Ebind 


tree 


*esc 


tool 


name 


(10 5 M ) 


(pc) 


(dex) 


(M @ ) 


(10 45 ergs) 


(10 3 yr) 


(10 5 yr) 


(10 s yr) 


47 Tuc 1 


15 


2.8 


-0.8 


0.1 


0.5 


3 


1-2 


0.5 - 1 


M2 2 


8.8 


0.9 


-1.6 


3 


100 


6 


0.4 - 4 


0.07 - 1 


M15 2 


5 


6.4 


-2.3 


0.3 


40 


10 


3 - 10 


1-5 


NGC 6624 3 


3.3 


0.8 


-0.4 


0.02 


0.5 


15 


0.3 - 0.4 


2 



iFreire et sd] j200l| l 

van Loon et al~] ]2006fr 
^Grindlay fc Liller| ^977\ 



effective. Though no challenge energetically, the 
100 pc dimension means that the escape timescale 
is so long (> 10 5 yrs) that much of the ICM would 
still be present when the next nova occurs. These 
systems are therefore much closer to the original 



limit discussed by Scott & Durisen ( 1978 ) , which 



predicts a wind leaving the system at the expected 
A^icm, but with a kinetic energy flux given by 
E n0 va and the CNe recurrence time. Such a wind 
would have a terminal velocity of « 60 km s _1 
(independent of the dSph mass), and a density at 
100 pc of - 10- 7 cm- 3 for a dSph with 1O 4 M of 
old stars. 
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